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We investigate Lagrangian relative dispersion in direct numerical simulation of two-dimensional 
inverse cascade turbulence. The analysis is performed by using both standard fixed time statistics 
and an exit time approach. Our results are in good agreement with the original Richardson's 
description in terms of diffusion equation. The deviations are only of the quantitative nature. 
CV These deviations, and the observed strong sensitivity to the finite size effects, are associated with 

the long-range-correlated nature of the particles' relative motion. The correlation, or persistence, 
^■* ' parameter is measured by means of a Lagrangian "turning point" statistics. 
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QJ ■ I. INTRODUCTION 

Understanding the statistics of particle pairs dispersion in turbulent velocity fields is of great interest for both 
theoretical and practical implications. At variance with single particle dispersion which depends mainly on the large 
scale, energy containing eddies, pair dispersion is driven (at least at intermediate times) by velocity fluctuations 
at scales comparable with the pair separation. Since these small scale fluctuations have universal characteristics, 
independent from the details of the large scale flow, relative dispersion in fully developed turbulence is expected 
to show universal behavior |l|,|2|. From an applicative point of view, a deep comprehension of relative dispersion 
mechanisms is of fundamental importance for a correct modelization of small scale diffusion and mixing properties. 

Since the pioneering work by Richardson |3(] , many efforts have been done to confirm experimentally or numerically 
his law [EpHl(J. Nevertheless, the main obstacle to a deep investigation of relative dispersion in turbulence remains 
the lack of sufficient statistics due to technical difficulties in laboratory experiments and to the moderate inertial 
range reached in direct numerical simulations. 

In this Paper we present a detailed investigation of the statistics of relative dispersion from extensive direct numerical 
simulations of particle pairs in two-dimensional Navier-Stokes turbulence. We will see that the main ingredient of 
the original Richardson description, i.e. Richardson diffusion equation, is sufficient for rough description of relative 
dispersion in this flow. Nevertheless, our simulations show that two-particle statistics is rather sensible to finite size 
J> . effects. This asks for a different analysis based on doubling time statistics which has recently introduced for the 
analysis of relative dispersion (llj. Comparison of numerical results with ones based on the Richardson's equation 
show that the last delivers a qualitatively good description of the doubling-time-distributions. The quantitative 
deviations found are attributed to the fact that the dispersion process is not purely diffusive and is influenced by 
ballistic (persistent) motion. 

The article is organized as follows: In Section || we discuss the Richardson's approach to the two-particle dispersion, 



> 



o 
o 
o 



c 



in Section III the fixed-scale properties of dispersion process (such as doubling-time statistics) are considered. The 
numerical approach and the results of simulations are discussed in Section [V. Section M is devoted to conclusions. 



The mathematical details of calculations of doubling-time statistics for the Richardson's case are given in Appendices 
and [|. 

II. STATISTICS OF RELATIVE DISPERSION 

Relative dispersion in turbulence is often phcnomcnologically described in terms of a diffusion equation for the 
probability density function of pair separation p(r, t) 



d P (r,t)_ d ( K{Tt) dpptr ] (1) 



dt drj \ ' dr,j 

with a space and time dependent diffusion coefficient K(r,t) |2J]. The original Richardson proposal obtained from 
experimental data in the atmosphere corresponds to K{r,t) = K{r) = k^e 1 ' 3 r 4 ' 3 B, where e has the dimension of 
energy dissipation (see below) and kg is an adimensional constant. In the d-dimensional isotropic case, this diffusion 
equation takes the form 

dp(T,t) _ 1 d d _ x dp{T,t) 

dt - r ^dr [) 8r () 

Its solution leads to the well known non-Gaussian distribution 

, ^ A ( 9rV3 \ 

*''*) = JWe^{-4k^m) (3) 

where A is a normalizing factor. The growth of pair separation is described in terms of a single exponent as 

R 2n (t) = (r 2n (t)) = C 2n e n t 3n . (4) 

The so called Richardson constant in (||) is g — C 2 = ^^k\. 

The Richardson's conjecture was formulated based on the scaling nature of the diffusion coefficient and in analogy 
with diffusion. No information about the nature of the pdf was available by that time. Any chose of the form 
K(r,t) ~ r 4 / 3 ~ Q (rf(£)) a / 2 (i.e. K(r,t) ~ r */3-<*t 3a / 2 ) wou i ( j give t } ie same sca ii ng i aw b 2 cx t 3 but with different 
pdf's (see @,||f2]jl|). 

The possibility to describe the dispersion process by means of a diffusion equation is based on essentially two 
important physical assumptions which can be verified a posteriori. The first one is that the dispersion process is 
self-similar in time, which is probably true in non intermittent velocity field M, the second one is that the velocity 
field is short correlated in time pj. In the limit of velocity field (5-correlated in time the diffusion equation (fil) 
becomes exact ]I5| , ^6| . As we proceed to show, the Richardson's conjecture (||), which is exact under small values 
of the persistence parameter of the flow, vide infra and [p.4| , |l7[ , still delivers a qualitatively good approximation for 
realistic 2d turbulent flows, whose persistence parameter of the order of I. 

Richardson scaling in turbulence is a consequence of Kolmogorov scaling for the velocity differences |2J. Under 
Kolmogorov scaling, the mean square relative velocity and the correlation time in the inertial range are given by 

(5v(rf) = vl (±\ ~ e 2 l 3 r 2 ' 3 (5) 



and 



r(r)=ro(^) 2/3 ^s- 1 / 3 r 2 / 3 ) (6) 

where rg, tq and vq are some (large scale) characteristic length, time and velocity scale and e ~ v 2 /tq is the energy 
flux in the inertial range. The value of the dimensionless combination Ps = ^ 3 ^ 2 - remains however unspecified by 
scaling considerations. It is referred to as a persistence parameter of the flow and plays a central role in describing 
single particle diffusion and pair separation | fi4] , p7t . 

The persistence parameter gives the ratio of the velocity correlation time to the Lagrangian characteristic time. In 
order to see how Ps influences Lagrangian dispersion, let us consider the following simple model, which has been used 
as a basis for building a stochastic model of turbulent dispersion |lj| . We take that the magnitude of the separation 
velocity is a function of r only so that Sv(r) = Vo(r /vq) 1 ' 3 . The temporal changes of the flow can be accounted for by 
letting the particle change its velocity direction from time to time, while keeping the velocity's magnitude constant. 
The probability that the relative velocity of particle separation changes its direction during time interval dt can be 
assumed to be dp — dt/r(r). The growth of the magnitude of the interparticle separation r(t) in dt is dr ~ 5v(r)dt, 
thus the probability to change the direction of velocity within dr is, using (g[|6|) 

dp = p(r)dr = ———— = — = — . (7) 

ov{r)T(r) VqTq r Ps r 

The distribution of the position of turning points in the separation follows from (m) [jf7| . The conditional probability 
density to find a next turning point at r 2 provided a previous one was at n < r 2 is given by 



*(r 2 |n)=p— (-) • (8) 

rsn V r i/ 

Note that the dependence of ^(r2\r\) only on the relative positions of the turning points, i.e. on t^/ti, is a clear 
consequence of scaling. 

The tail of ^(r2\r\) decides about the existence of the second moment of this distribution, i.e. on the fact whether 
the corresponding motion is short- or long-range correlated in space. Depending on the persistence parameter Ps, 
the dispersion can be either diffusive (Ps <C 1) or ballistic (Ps ^ 1) in nature. Note that the power-law tail of the 
distributions make the problem extremely sensitive to the finite-size effects, especially for large Ps, when the weights 
of ballistic events (Levy- walks Jig] ) is considerable. 

We note that the value of Ps is not a free parameter, but is fixed for a given physical situation. The scaling nature 
of turbulence supposes that this parameters is a constant, depending only on general properties of the flow, e.g. on 
its 2d or 3d nature (in this last case also the overall geometry of the flow can be of importance). On the other hand, 
since the nature of dispersion process depends crucially on the value of Ps, the only way for getting quantitative 
information about the dispersion is through direct numerical simulations. The strong finite-size effects in relative 
dispersion simulations require the introduction of quantities which are less sensitive to finite resolution. 

III. EXIT TIME STATISTICS 

All the results discussed in the previous Section can be observed only in high Reynolds number flows in which the 
inertial range, where the scaling laws (|5|,|6|) hold, is sufficiently wide. The needs for large Reynolds is particularly 
severe in the case of Richardson dispersion, as a consequence of the long tails in (0). Moreover, the observation of 
time scaling laws as (Q) requires sufficiently long times in order to forget the initial separation r(0) U. 

For these reasons, the observation of Richardson scaling (i.e. (0) or (0)) is very difficult in direct numerical 
simulations where the Reynolds number is limited by the resolution. The same kind of limitations arise in laboratory 
experiments, as a consequence of the necessity to follow the Lagrangian trajectories which limits again the Reynolds 
number |0,||. 

To partially overcome these difficulties, an alternative approach based on doubling time (or exit time) statistics has 
been recently proposed P,|fl|. Given a set of thresholds R n — p n R(0) within the inertial range, one computes the 
"doubling time" T p (R n ) defined as the time it takes for the particle pair separation to grow from threshold R n to the 
next one R n +i- Averages are then performed over many dispersion experiments, i.e., particle pairs, to get the mean 
doubling time (T p (R)). The outstanding advantage of this kind of averaging at fixed scale separation, as opposite to 
a fixed time, is that it removes crossover effects since all sampled particle pairs belong to the same scales. 

The problem of doubling time statistics is a first-passage problem for the corresponding transport process. For 
the Richardson case, in 2d it is given by the solution of the Richardson's diffusion equation, Eqs. (0), with initial 
condition p(r, 0) = 6(r — R/p)/2n and absorbing boundary at r = R (so that p(R, t) = 0). The pdf of doubling time 
can be obtained as the time derivative of the probability that the particle is still within the threshold 
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P°^ = -T+ P(r,t)dr. (9) 

at J\r\<R 



Using (g) one obtains 

(10) 



p D (t) = -2^/3^4/3 ^M 



r=R 



The solution using the eigenfunction decomposition is given in Appendix |A| and shows that the long-time asymptotic 
of pr> (t) is exponential: 



PD 



(t) ~ exp(-Kfc £ 1/3 i?~ 2/3 £) (11) 



where k ss 2.93 is a number factor. This exponential nature of the tail of p£>(£)-distribution will be confirmed by 



direct simulations in Section IV 



Note that the combination e _1 / 3 i? 2 / 3 has a dimension of time and is proportional to the average doubling time 
(Tp(R)}. This time can be obtained by a simple argument reported in the Appendix |B|. In the 2-dimensional case one 
obtains 



3 p 2 ' 3 - 1 R 2 ' 3 

< T ^ = 4&7v7^— (12) 

Prediction (n2[) contains the parameter fc which is, as shown in Section O, is equivalent to the Richardson constant 
g. As a consequence, the computation of average doubling time can be used for an alternative (and more robust, as 
we will see) estimation of the Richardson constant. It is convenient to rewrite the doubling time pdf ([fl]) in terms of 
the average doubling time (T p (R)}. Making use of ( |l2|) one obtains in 2d 

PD (t) ~ exp(-0.252 t/(T(R)}) . (13) 

which is a parameterless, universal function. 

IV. DIRECT NUMERICAL SIMULATIONS 

Pair dispersion statistics has been investigated by high-resolution direct numerical simulations of the inverse energy 
cascade in two-dimensional turbulence |L9 |. There are several reasons for considering 2D turbulence. First of all, 
the dimensionality of the problem makes feasible direct numerical simulations at high Reynolds numbers. Moreover, 
the observed absence of intermittency [£0[ makes the 2D inverse energy cascade an ideal framework for the study of 
Richardson scaling in Kolmogorov turbulence. 

The 2D Navier-Stokes equation for the vorticity w = Vxv = —Aip is: 

dtuj + J (w, ip) — isAui — atu + cf), (14) 

where tp is the stream function and J denotes the Jacobian. The friction linear term —auj extracts energy from the 
system to avoid Bose-Einstein condensation at the gravest modes [^IJ. The forcing is active only on a typical small 
scale // and is (5-correlated in time to ensure the control of the energy injection rate. The viscous term has the role 
of removing enstrophy at scales smaller than If and, as customary, it is numerically more convenient to substitute it 
by a hyperviscous term (of order eight in our simulations). Numerical integration of (U4) is performed by a standard 
pseudospectral method on a doubly periodic square domain at resolution TV = 1024 or TV = 2048. All the results 
presented are obtained in conditions of stationary turbulence. 

In Figure p] we plot the energy spectrum, which displays Kolmogorov scaling E(k) = Ce 2 / 3 k~ b l 3 over about two 
decades with Kolmogorov constant C ~ 6.0. In the inset we plot the third order structure function S3 (r) — (Sv(r) 3 ) 
compensated with the theoretical prediction S3 (r) = 3/2er. The observation of the plateau confirms the existence 
of an inverse energy cascade and indicates the extension of the inertial range. Previous numerical investigation has 
shown that velocity differences statistics in the inverse cascade is not affected by intermittency corrections [EOf. In 
this case we may expect the Lagrangian statistics to be self-similar with Richardson scaling Kj . 

Lagrangian statistics is obtained by integrating the trajectories of 64000 particle pairs in the turbulent velocity 
field, initially uniformly distributed with constant separation R(0). 

The Lagrangian data reported below are in dimensionless units in which separations are rescaled with the box size 
L = 2-7T and time with the large scale time To = (L 2 /e) 1 / 3 . 

A. Relative dispersion analysis 

In Figure || we plot the relative dispersion R 2 (t) for two different initial separation, R(0) = 5x/2 and R(0) = Sx 
(where Sx = 2n/N is the grid mesh). The Richardson t 3 law (H) is observed in a limited time interval, especially for 
the larger i?(0) run. It is remarkable that the relative separation law displays such a strong dependence on the initial 
conditions even in our high resolution runs. 

This dependence makes the determination of the Richardson constant particularly difficult. In Figure we also 
show the compensated plot R 2 (t)/t 3 which, in the adimensional units, should directly give the constant g. It is clear 
that a precise determination of g is impossible; even the Richardson scaling (H), when looked in a compensated plot, 
is rather poor. Figure |2j suggests that starting with an intermediate initial separation would give a wider scaling 
range. Of course, one would like to avoid this "fine tuning" , which is probably impossible to implement in the case of 
experimental data. 

The probability distribution function of pair separations is plotted in Figure for the R(0) — 5x12 run. At short 
time t — 0.015, in the beginning of the t 3 range in Figure g, we found that the Richardson pdf (||) fits pretty well 
our data, although some deviations can be detected. Of course, at time comparable with the integral time t = 0.77, 



particle separations are of the order of the integral scale and we observe Gaussian distribution. The crossover between 
these two regimes is extremely broad: deviations from Richardson pdf are clearly seen already for the times well 
within the Richardson's t 3 range. To observe better this transition, in Figure |I] we plot, in log-log plot, the right tail 
of — ln(p(r, t)/p(0, £)). We see that the tail slope can be fitted with an exponent a which change continuously in time, 
from 2/3 to the Gaussian value 2 (see the inset of Figure 0). Thus self similarity, if it exists, is reduced to the very 
short time at the beginning of dispersion. Moreover, the scaling region is strongly affected by the choice of initial 
separation, as shown in Figure 0. 

B. Doubling time data 

From the same Lagrangian trajectories discussed in the previous Section, we have computed the doubling time 
statistics. The average doubling time as function of the scale is plotted in Figure [|. The line represent the dimensional 

prediction (T(R)) ~ i? 2 / 3 . In the inset we plot the quantity ^ — — r-*— e £L 9 which, from ([|) and ( |l2| ) gives the value 
of the Richardson constant g ~ 3.8, which is compatible with the estimation from Figure 0. It is interesting to compare 
our result with previous estimations of g. The only experimental estimation of g for 2D inverse cascade m gives a 
value about 7 times smaller (but the Reynolds number in the experiment is even smaller than in present simulations). 
Other estimations of g are based on kinematic simulations with synthetic flows. In all these cases J5plJ22| the reported 
values are even smaller. On the other hand, our numerical finding is not very far from the prediction of turbulence 
closure theory |1^J23]] . 

From Figure g we observe that at very small separations R ~ 10 -3 , the doubling time has a tendency to a 
constant value {T(R)} ~ 0.0016. On these scales we are below the forcing scale (see Figurey]), and the velocity 
field can be assumed smooth. As a consequence of Lagrangian chaos we expect on these scales an exponential 
amplification of separations |J24| at a rate given by the Lagrangian Lyapunov exponent A. The latter can be obtained 
as A = linifl^o hi p/(T(i?)) Wu and gives A ~ 110 (in adimensional units). The Lagrangian Lyapunov exponent A is 
a small scale quantity (i.e. depends on the Reynolds number of the simulation), and thus has to be compared with 
a small scale characteristic time. One can estimate the smallest characteristic time T m i n by the minimum value of 
{k^Eik))- 1 . We obtain A ~ 0.23r"^. 

The comparison of Figures and pf shows the advantages of fixed scale statistics with respect to fixed time statistics. 
In particular, the dependence of R(t) and the Richardson constant on the initial separation are absent in fixed scale 
statistics T(R). 

In Figure ^| we plot the doubling time pdf pd(T) compensated with the mean value {T(R)) at different scales in 
the inertial range 0.003 < R < 0.046. First, we obtain a very nice collapse of the different curves, indicating that the 
process is really self-similar. Second, we observe the exponential tail predicted in Section fill with a fitted coefficient 0.3 
which is indeed not far from the theoretical prediction (|l|) based on the Richardson's picture. The difference between 
the predicted and measured values of the prefactors is not large, but perceptible: it shows that the Richardson's 
equation gives a correct qualitative description of the dispersion process, but is not exact. The reasons for deviations 
from the diffusive picture proposed by Richardson are the long-range correlations in the particles' motion, as seen 
from the analysis of the turning points of their relative trajectories. 

C. Turning points statistics and the persistence parameter 

A possible explanation for the deviations of our high resolution numerical data from the Richardson's picture is the 
not too small value of the persistence parameter. As discussed in Section O, at large values of Ps the contribution 
of ballistic events may lead to non-Richardson distributions and moreover makes the dispersion strongly sensible to 
finite-size effects, cutting the longer trajectories. 

We have computed the persistence parameter making use of (g). We have recorded, for each pair, the set of turning 
points Ti at which the pair's relative velocity changes sign. From the set of r^ we have then computed the pdf of 
the ratio r i+ i/ri, accumulating for all the i and all the pairs. The result, plotted in Figure |?], gives Ps ~ 0.87. The 
requirement that both t\ and ri are in the inertial range, strongly limits the statistics on turning points and the 
numerical result is affected by rather large uncertainty. Nevertheless, it is remarkable that the power law tail in the 
conditional probability density ^(r^ri) is well observed in our numerical simulations. This justifies, a posteriori, the 
use of models based on 4 r (f2|ri) for describing relative dispersion 111]. The numerical value of the effective persistence 
parameter Ps ~ 0.87 is not so small, and can explain the observed deviations from Richardson pdf (which are however 



less pronounced in a 2d flow than in a theoretical one-dimensional model |17|): The transport in a 2d turbulent flow 
is neither purely diffusive nor ballistic. 



In order to be more confident on the numerical value of Ps obtained through the turning-points statistics, let us 
show that it agrees with a simple estimates based on the values of the Kolmogorov's and the Richardson's constants. 
According to the Kolmogorov's scaling, the mean squared relative velocity of the pair is given by 

(5v 2 (r))=C 2 e 2 / 3 r 2 /\ (15) 

with C2 — 12.9 p3]. If the particles separate ballistically with the rms velocity 

Sv(r) = C£/V/V/ 3 (16) 



the distance between them should grow as 



3 



i?L x = ( I ) cTet 3 • (17) 



On the other hand, due to the unsteadiness of the separation velocity, the distance between the particles grows slower, 
namely as R 2 — get 3 , so that the factor 



t-2 r>2 I r>2 I "\ 9 



i=RVK^=[^j -^ (is) 

serves as a measure of this unsteadiness and £ 2 is connected with the value of the persistence parameter. In our 
case, £ 2 w 0.28. Within the stochastic model of Ref. |lj| this corresponds to a value of Ps between 1.1 and 1.2, 
in reasonable agreement with the direct measurement from the turning-point statistics, and again corroborates the 
stochastic approach. 

We also note a possibility to "tune" the Ps value by performing simulations in which the Lagrangian trajectories are 
integrated according to x = Av(x, t). By changing the value of parameter A one effectively changes vq and thus Ps. 
In the extreme case A — > the trajectories resemble those in a time ^-correlated velocity held. In the opposite limit, 
A >> 1 we have dispersion on a quenched field. Of course, it is only for the standard value A = 1 that Lagrangian 
trajectories move consistently with velocity field (i.e. for v = a = (\> = (|14|) conserves vorticity along the Lagrangian 
trajectories). For other values of A such simulations suffer the typical problem of advection in synthetic field (i.e. 
wrong reproduction of the sweeping effect, see |9j for a discussion). Simulations for several values of A show that 
existence of the power- law tails of ^(^Irj) is a robust effect, as supposed by the model of Refs. p^Jl7| , and that Ps 
grows with A. As an example, in Figure |we plot the probability density ^(^Iri) obtained from a simulation with 
A = 0.5. All the Eulerian parameters are the same of Figure ul We again observe a clear power law tail but now with 
Ps~0.58. 

V. CONCLUSIONS 

In this paper we have investigated the Lagrangian relative dispersion in direct numerical simulation of two- 
dimensional turbulence. The inverse energy cascade of two-dimensional turbulence displays Kolmogorov scaling 
without intermittency and it is thus the natural framework for investigating possible deviations from the classical 
Richardson picture. Moreover, the large enough values of the accessible Reynolds numbers make it possible to obtain 
the quantitative results not inferior to those of laboratory experiments. 

The analysis of the numerical data was performed by using both standard statistics at fixed time and exit time 
statistics at fixed scale. The latter is shown to be more robust to finite Reynolds effects. An application of exit time 
statistics is developed for measuring the Richardson constant. The value found is somewhat larger than one obtained 
experimentally and in kinematic simulations, and thus nearer to the predictions of the closure theories. 

The results of numerical simulation of statistics of the intcrparticle distances are in good agreement with the 
original Richardson's description in terms of diffusion equation. The large deviations (with respect to Richardson 
theory) observed in the pdf of separation at fixed times are mostly related to crossover effects due to finite Reynolds 
numbers and disappear when looking at exit time statistics. The deviations found here are only of the quantitative 
nature. Thus, the Richardson's equation gives a good basis for qualitative description of the dispersion in turbulent 
flows. 

The deviations found in the fixed-scale statistics and also the strong sensitivity to the finite size effects can be 
associated with the long-range-correlated nature of the particles' relative motion in turbulent flow. Paying attention 
to the turning points of the relative trajectories allows for estimating the effective persistence parameter of the motion 



which is found to be of the order of unity. Thus, the motion shows a strong ballistic component and is not purely 
diffusive. On the other hand, the correlations are not too strong to fully destroy the Richardson's picture. 

We note that the methodology of analysis proposed here based on the fixed-scale statistics and on the analysis of the 
relative trajectories can be also applied to the analysis of laboratory experiments. It would be extremely interesting 
to see, to what extent the laboratory flows reproduce the results reported here. 
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APPENDIX A: THE PDF OF DOUBLING TIMES 

Let us discuss the probability pu (t) of the time when the a pair of particles initially at distance R/p separates up 
to the distance -R, and obtain its asymptotic decay of this probability for t large. 

Changing to a variable £ = (koe 1 ^ 3 )^ 1 ^ 2 ^^ 3 reduces (||) to a radial part of a spherically symmetric diffusion equation 
in 3d dimensions with constant diffusion coefficient. In 2d one has 

^ - — — £ 5 —v (AD 

with the initial condition p(£,0) = <$(£ m in — £) with £ m ; n = {kQS 1 ^ 3 )^ 1 ^ 2 {R/ p) 1 ^ 3 and with the boundary condition 

P(£max, t) = 0, With £ max = (fc £ 1 /3)-l/2 i? l/3. 



The solution of a boundary- value problem for (Al) can be obtained by means of eigenfunction decomposition. As- 
suming the variable separation we get the solution in the form p(£, t) = ^\ e _Ai *'0i(£)j where t/^(£) is an eigenfunction 
of the equation 

satisfying the boundary condition -0(£ max ) = 0. The corresponding solution which is nonsingular in zero is ipi = 
£~ 2 J2(3Ai£) (J2 is the Bessel function p5j). The fact that ip vanishes at £ max gives 3\i£, m ax = J24, where ji,i is 
the «-the real zero of J2(x). For example, the smallest eigenvalue is X 2 = j| i/9£max ~ '2-Q'&koe 1 / z R~ 2 / 3 . Thus, the 
long-time asymptotic of the doubling-time distribution is exp(— 2.93koS 1 / 3 R~ 2 / 3 t). 

APPENDIX B: AVERAGE DOUBLING TIME 

The mean doubling time can be obtained from a stationary solution of the Richardson diffusion equation. Imagine 
that one particle per unit time is introduced at r = R/p and there are respectively a reflecting and absorbing 
boundaries at r = and r = R. The stationary solution of (|2]) in 2d with the appropriate boundary conditions and 
continuity are R/p is 



p(r) 
The number of particle in r < R is 



3 4 / 3 - I] for < r < R/p 



-4/3 



£) '°-l for R/p<r<R 



(Bl) 



N = I p(r)dr = 2ir / rp(r)dr (B2) 

r\<R 



where Sd is a surface of a unit sphere in d dimensions. By using (Bl) one obtains 



N = 2nC(l - p- 2/3 )R d . (B3) 



The current at r = R, i.e. the number of particle exiting from the boundary R per unit time, is given, as in (QOh , as 

8tt 



J = -2^/3 fcoi? 7/3 9p(r) 



r=R 3 



Ce 1/3 k R^ 3 (B4) 



Or 
The mean doubling time is the average time spent by a particle at r < R. It is given by the ratio N/J and thus 



(T p (R)) = l f/ 3 lj /3 R 2 / 3 (B5) 



which is Eq.ji"^). 
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FIG. 1. Energy spectrum of the inverse cascade simulations E(k). The dashed line has the Kolmogorov slope —5/3. In the 
inset it is shown the compensated third order longitudinal structure function Ss(r)/(er). The dashed line is the prediction 3/2. 
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FIG. 2. Relative dispersion R 2 (t) with R(0) — 5x/2 (+) and R(0) — 8x (x). In the inset the compensated plot R 2 (t)/(et 3 
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FIG. 3. Probability distribution function of relative separations at times t — 0.031 (*) and t = 0.77 (x) rescaled with 
R(t) = (r 2 (i)) 1 ' 2 . The continuous line is the Richardson prediction (pi), the dashed line is the Gaussian distribution. 
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FIG. 4. Right tail of -log(p(r,t)/p(0,t)) at times t = 0.015 (+), t = 0.041 (x), t = 0.067 (*) and t = 0.77 (□) in log-log 
plot. The two lines represent the Richardson slope 2/3 and the Gaussian slope 2. The inset shows the exponent of the right 
tail of the pdf as a function of time. 
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FIG. 5. Mean doubling time (T(R)} as function of the separation R. The ratio is r = 1.2 and the average is obtained over 
about 5 x 10 5 events. The line represent the dimensional scaling R 2 ' 3 . In the inset the compensated plot gives the value of 
Richardson constant g ~ 3.8, as explained in the text. 
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FIG. 6. Pdf of doubling times for distances R = 0.003 (+), R = 0.075 (x), R = 0.02 (*) and R = 0.046 (□). The dashed 
line is the exponential exp(— 0.3 T/{T)). 
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FIG. 7. Probability density function of turning point ratio ^/(ra/ri). The exponent of the power law (dashed line) gives 
the value Ps ~ 0.87. 
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FIG. 8. The same of Figure M but for A = 0.5. The persistence parameter is now Ps ~ 0.58. 
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